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ABSTRACT 



Aims. We present modeled flux and linear polarisation signals of starlight that is reflected by spatially unresolved, 
horizontally inhomogeneous planets and discuss the effects of including horizontal inhomogeneities on the flux and 
polarisation signals of Earth-like exoplanets. 

Methods. Our code is based on an efficient adding-doubling algorithm, which fully includes multiple scattering by gases 
and aerosol/cloud particles. We divide a model planet into pixels that are small enough for the local properties of the 
atmosphere and surface (if present) to be horizontally homogeneous. Given a planetary phase angle, we sum up the 
reflected total and linearly polarised fluxes across the illuminated and visible part of the planetary disk, taking care to 
properly rotate the polarized flux vectors towards the same reference plane. 

Results. We compared flux and polarisation signals of simple horizontally inhomogeneous model planets against results 
of the weighted sum approximation, in which signals of horizontally homogeneous planets are combined. Apart from 
cases in which the planet has only a minor inhomogeneity, the signals differ significantly. In particular, the shape of the 
polarisation phase function appears to be sensitive to the horizontal inhomogeneities. The same holds true for Earth-like 
model planets with patchy clouds above an ocean and a sandy continent. Our simulations clearly show that horizontal 
inhomogeneities leave different traces in flux and polarisation signals. Combining flux with polarisation measurements 
would help retrieving the atmospheric and surface patterns on a planet. 

Key words, exoplanets - polarisation - numerical code 



1. Introduction 



Since [Mayor fc Quelozl (|1995f ) discovered the first planet or- 
biting another main sequence star almost two decades ago, 
the rapid improvement of detection methods and instru 
merits has yielded hundreds of exoplanets, in cluding sev 
eral te ns of so-called sup e r-Ear t hs (see e.g. iLeger et al 



2011 



Charbonneau et all l2009t iMiller-Ricci fc Fortnev 



2010; IBeaulieu et all 12006). and many more will follow in 



the coming years. The next step in exoplanet research is 
the characterisation of the atmospheres and surfaces (if 
present) of detected exoplanets: what is their composition 
and structure? 

Currently, exoplanet atmospheres are being charac- 
terised using the transit method (see e.g. IBeaulieu et al.l 
l2010t IMiller-Ricci fc Fortnevll2010D . This method is based 
on measurements of the wavelength dependendence of 
starlight that filters through the upper planetary atmo- 
sphere during the primary transit, or of the planetary 
flux just before or after the secondary eclipse. The tran- 
sit method is mostly applied to gaseous planets that orbit 
close to their star. The chances to catch gaseous planets 
in wide orbits, such as Jupiter and Saturn in our Solar 
System, transiting their star are extremely small, because 
their orbital plane should be perfectly aligned with our line 
of sight and because their transits are very rare. Earth- 
sized exoplanets in the habitable zone of a solar-type star 
are probably too small and transit too seldom to reach 



a sufficient signal-to-noise r atio to do transit spectroscopy 
(jKaltenegger fc Traubll2009f) . 

A promising method to characterise atmospheres and 
surfaces (if present) of exoplanets that are small and/or in 
wide orbits, is direct detection, in which the starlight that a 
planet reflects and/or the thermal radiation that a planet 
emits is measured separately from the stellar light (except 
for some background starlight). Some instruments that are 
being designed for such direct detections are SPHERE (for 
the VLT) and EPICS (for the European Extremely Large 
Telescope, or E-ELT). An example for a space telescope 
for direct detection is the New Worlds Observer (NWO) 
(jCash k New Worlds Study Teaml l2010h . which is under 
study by NASA. Through direct detections, broadband im- 
ages and/or spectra of various types of exoplanets will be- 
come available in the near future. 

Knowing the Solar System planets, it is to be expected 
that exoplanets that will be observed are horizontally in- 
homogeneous, e.g. with clouds and hazes in patches, such 
as on Earth and Mars, or in banded structures, such as, for 
example, on Jupiter and Saturn. And there will undoubt- 
edly be solid exoplanets with significant local variations in 
surface reflection properties and texture, such as the Earth 
with its continents and oceans. Although the lack of spatial 
resolution in near future exoplanet observations will merge 
all spatial variation into a single image pixel and/or spec- 
trum, accounting for the existence of horizontal inhomo- 
geneities will be important when trying to retrieve planet 
characteristics. For example, efforts to identify spectral sig- 
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natures of life on other planets will face various challenges, 
such as clouds masking or m imicking the signatures of veg- 
etation (|Tinetti et al.l 1200(1 and references therein). The 
cloud coverage will also influence the retrieval of mixing 
ratios of atmospheric gases, such as water vapour and oxy- 
gen, from reflected light spectra. In particular, the larger 
the fraction of clouds across an exoplanetery surface, the 
smaller the depth of gaseous absorption bands. Absorption 
band depths are, however, also influenced by the cloud top 
altitudes, with higher clo uds yielding shallower absorption 
bands (for examples, see IStaml |2008|) . The distribution of 
cloud top altitudes across a planet will thus also be a pa- 
rameter to take into account. 

Horizontal inhomogenities can have large effects on the 
flux of starligh t that is reflected by a planet. In particu- 
lar, iFord et all (|200lf ) have shown that an Earth analogue 
planet without clouds would show diurnal flux variations of 
up to 150% due to the variation of the al bedo of different 
region s on the planet. In full agreement, lOaklev fc Cash! 
(2009) have calculated that in absence of an atmosphere, 
the flux of an Earth-like exoplanet would show a clear di- 
urnal variability as different continents would rotate in and 
out of the field of view of an observer. The presence of 
clouds in the Earth's atmosphere significantly complicates 
the ch aracterisation of various surface types. IFord et all 
(|2001|) show, for example, that an Earth-like cloud pat- 
tern would s uppress the diurnal flux variations to as little 
as 20%, while lOaklev fe Ca sh (2009) indicates that charac- 
terizing the surface of the planet is possible only for cloud 
coverages significantly lower than the average coverage on 
Earth (< 25% versus ~ 60%). To quantitatively estimate 
the effects of horizontal inhomogeneities due to clouds or 
surface features on observed spectra, and to be able to ac- 
count for such variations in the retrieval of planet char- 
acteristics from future observations, numerical codes are 
essential tools. 

In this paper, we present our numerical code to calculate 
spectra of starlight that is reflected by spatially unresolved, 
horizontally inhomogeneous exoplanets. The main differ- 
ence with other codes for horizontally inhomogeneous plan- 
ets (s uch as those u sed bylFord et al.ll200lt[Oaklev fc Cash! 
2009| iTinetti et all 120061 ) is that it can be used to calcu- 
late not only the flux of reflected starlight but also its 
state (degree and direction) of polarisation. Polarimetry 
promises to play an important role in exoplanet research 
both for exoplanet detection and characterisation. In par- 
ticular, because the direct starlight is unpolarized, while 
the starlight that is reflected by a planet will usually be 
polarized Csee e. g. IZugger et al.l 201C Stam et all l2006t 
120041: ISaar fc Seaeerll2003l: ISeager et aUhoOOft . polarimetry 
can increase the pl anet-to-star contr ast ratio by 3 to 4 or- 
ders of magnitude (|Keller et al.ll2010l) , thus facilitating the 
detection of an exoplanet that might otherwise be lost in 
the glare of its parent star. Polarimetry will not only help to 
detect a planet, it will also confirm the status of the object, 
since background objects will usually be unpolarized. 

The importance of polarimetry for studying plane- 
tary atmospheres and surfaces has been shown many 
times using observations of th e Earth and other Solar 
System planets (see for example Hansen fc Hovenier 19741 
Hansen fc Travis! Il974t iMishchenkd Il990t iTomasko et ajj 



ISeaeer et al.l l2000h . In particular the sensitivity of polar- 
isation to the microphysical properties of the scatterers in 
the planetary atmosphere, make it a crucial tool for braking 
degeneracies that flux only observations can have. 

The radiative transfer calculations in our code are based 
on an efficient adding-doubling algorithm (|de Haan et al.l 
11987! ) which fully includes multiple scattering by gases and 
aerosol/cloud particles, that was used before for flux and 
polarisatio n calc u lation s for gaseous and t errestr i al exo - 
planets bv IStaml (1200311 : IStam et al.l (12004ft : IStaml (|2008h : 
iKaralidi et al.l ()2011 ). These authors, however, assumed 
each exoplanet to be horizontally homogeneous, such that 
it could be treated as a single starlight scattering 'particle', 
which allowed for a very fast integration of the reflected 
flux and polarisation signals across the pl anet's disk for the 
whole planetary phase angle range (see St am et al.ll2006l 
for a description of this disk- integration algorithm). With 
this horizontally homogeneous code, the signals of horizon- 
tally inhomogeneous planets can be simulated using the so- 
called weighted sum approximation: signals of homogeneous 
planets are multiplied by a weighting factor and summed 
to yield the final signal. We, on the other hand, divide a 
horizontally inhomogeneous model planet into pixels that 
are small enough for the local properties of the atmosphere 
and surface (if present) to be horizontally homogeneous. 
For each type of pixel, we perform adding-dou bling radia- 
tive transfer calculations ()de Haan et al.l fl987h and, given 
the planetary phase angle, we sum up the reflected total 
and polarised fluxes across the illuminated and visible part 
of the planetary disk. Our code for horizontally inhomo- 
geneous planets allows investiging the applicability of the 
weighted sum approximation, and the effects of horizontal 
inhomogeneities on the flux and in particular the polarisa- 
tion signals of exoplanets. 

This paper is organised as follows. In Sec.O we describe 
our numerical method to calculate the flux and polarisation 
of starlight that is reflected by a horizontally inhomoge- 
neous planet, and in Sec. [31 we present simulations of flux 
and polarisation for different types of horizontal inhomo- 
geneities. In Sec.|4l we present flux and polarisation signals 
of horizontally inhomogeneous Earth-like planets and com- 
pare them to signals for horizontally homogeneous planets. 
Finally, in Sec. [5l we discuss and summarize our results. 
Appendix A contains the results of testing our code for hori- 
zontally inhomogeneous planets against an existing code for 
horizontally homogeneous planets. 

2. Calculating reflected starlight 

Light can fully be described by a flux vector ttF, as follows 



ttF 



F 

Q 

u 

V 



(1) 



2009), as well as by mod eling of Solar System planets or 



giant exoplanets ( e. g. | Ma dhusudhan fc Burrows 



IStaml I2008L 120031: IStam et al.l 120041: ISaar fc Seager 



2012; 



2003; 



with ttF the total, ttQ and nil t he linearly and ttV the 
circularly polar ised fluxes (see e .g. lHansen fc Travisl 11974 
iHovenier et"aT1l2004J: IStamll2008l ). Parameters ttF, ttQ, ttU 
and irV depend on the wavelength A, and have dimensions 
W m _2 m _1 . Parameters irQ and ttU are defined with re- 
spect to a reference plane, for which we choose the plane- 
tary scattering plane, i.e. the plane through the centers of 
the planet, star and observer. Note that this plane is usu- 
ally not the same as the planetary orbital plane; only for 
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orbits that are seen edge-on, the two planes coincide at all 
phase angles. In the following, we will ignore irV, because it 
is usually very small ([Hansen fc Travislll974l) . and because 
the errors in calculated val ues of nF, 7rQ, and nil due to 
ignoring ttV are neglig ible IjStam fc Hovenierl[2005h . 

The degree of polarisation P of vector ttF is defined as 



P 



y/Q 2 + U 2 
F 



(2) 



which is independent of the choice of reference plane. For 
planets that are mirror-symmetric with respect to the plan- 
etary scattering plane (i.e. horizontally homogeneous plan- 
ets), Stokes parameter U equals zero. In that case, we can 
use an alternative definition of the degree of polarisation 



P = — — 



(3) 



with the sign indicating the direction of the polarisation, 
i.e. if P s > (P s < 0) the light is polarized perpendicular 
(parallel) to the planetary scattering plane. The absolute 
value of P s is just equal to P. 

We calculate the flux vector of starlight that has been 
reflected by a spherical planet with radi us r at a distance d 
from the observer using (d > r) as (see lStam et~ai1l2006| ) 



nF(a) 



1 r 2 

--^S(a)irFo, 



(4) 



with a the planetary phase angle, i.e. the angle between 
the star and the observer as seen from the planet's center. 
Furthermore, ttFq is the flux vector of the incident starlight 
and S the 4x4 planetary scattering matrix. In the follow- 
ing, we normalise Eq. [H assumi n g r = 1 and d = 1. 

As described in IStam et aD ( 2006f l , the planetary scat- 
tering matrix S can be calculated by integrating local re- 
flection matrices R across the illuminated and visible part 
of the planetary disk, as follows 

S(a) = ~ [ MMoi(/3 2 )i2(M,Mo ) A0)i(/3 1 )dO ! (5) 

J D 

where dO is a surface element on the planet, and R is 
the local reflection matrix, which describes how starlight 
that is incident on a given location of the planet is re- 
flected towards the observer. The reference planes for R 
are the local meridian planes, which contain the direc- 
tion of propagation of the incident and reflected light, re- 
spectively, and the local vertical dir ection. The matrices 
L are so-called rotation matri ces (see IHovenier et al.ll20f)4 
IHovenier fc van der Medfl98l that are used to rotate from 
the planetary scattering plane to the local meridian planes 
and back: 



10 

cos 2/3 sin 2/3 

-sin 2/3 cos 2/3 

1 



(6) 



Angle /3 is measured in the anti-clockwise direction from the 
old to the new reference plane when looking in the direction 
of propagation of the light. 

In Eq. [5J we assume that the planetary atmosphere and 
surface (if present) are locally plane-parallel and rotation- 
ally symmetric with respect to the local vertical direction. 
Therefore, each matrix R depends on fi$ — cos 9$, with 8q 



the angle between the local zenith and the direction to- 
wards the star, on /i = cos 9, with 9 the angle between the 
local zenith and the direction towards the observer, and on 
A(f) = 4> — (j}Q, the angle between the azimuthal angles of the 
incident and the reflected light, respectively. Azimuthal an- 
gles are measured rotating clockwise when looking up from 
an arbitrary, local vertical plane towards the local vertical 
plane containing the direction of propagation of the light. 

Fo r horizontally homogeneous planets, IStam et ail 
(|2006l ) present an efficient method for evaluating Eq. [5] 
that uses an adding- doubling radiative transfer algorithm 
(|de Haan et alj|l987[ ) to calculate the coefficients of the ex- 
pansion of the local reflection matrix (which is the same 
across the planet) into a Fourier series. These coefficients 
are then used to compute coefficients of the expansion of 
matrix S into generalized spherical functions. With these 
expansion coefficients, S can be calculated rapidly for any 
phase angle a. Since with this method, a planet is basically 
treated as a single light-scattering particle, it cannot be 
used for horizontally inhomogeneous planets. 

To calculate S for horizontally inhomogeneous planets, 
we divide a planet into pixels small enough for the local 
atmosphere and surface (if present) to be considered both 
plane-parallel and horizontally homogeneous. For each type 
of pixel (a combination of surface and atmosphere proper- 
ties), we first calculate the coefficients of the expansion of 
the local reflection matrix R into a Fourier series , using the 
adding-doubling algorithm (jde Haan et al.ll987l) . Then, for 
each given planetary phase angle and each pixel, we use the 
respective Fourier coefficients to calcu late the local reflec- 
tion matrices (see Ide Haan et al.l[T987l) . The local matrices 
are summed up according to 

4 N 

S(a) = - S^ jiiHQiL(f} 2 i)Ri{ni,iioi, A(f>i)L(l3u)dOi, (7) 



with N the total number of pixels on the illuminated and 
visible part of the planetary disk. 

In the following, we assume unpolarised incident 
starlight, since integrated over the stellar disk, l ight of solar- 
type stars can be assumed to be unpolarised (|Kemp et al.l 
1987). In this case, ttFq = 7rFol, with 1 the unit column 
vector and ttFq the total incident stellar flux (measured 
perpendicular to the propagation direction of the light), 
for which we assume a normalized value of 1 W m~ 2 m . 
Thanks to the assumption of unpolarized incident starlight, 
rotation matrix L(f3n) can be ignored in Eq. [7] 

Because of the normalizations and the assumption of 
unpolarized incident light, the total flux that is reflected 
by a planet (cf. Eq. \£§ is given by 



nF n (X,a) = ian(A, a), 



(8) 



with an the (1,1) -e lement of matrix S (see IStaml [2008; 
iKaralidi et al.l [201 lh . The subscript n indicates the nor- 
malization. When a = 0°, the hence normalized total flux 
equals the planet's geometric albedo Aq. Our normalized 
total and polarized fluxes irF n , irQ n and TrU n can straight- 
forwardly be scaled to any given planetary system using 
Eq.[4]and inserting the appropriate values for r, d and ttFq. 
The degree of polarization (P or P s ) is independent of r, d 
and 7rPoi and would thus not require any scaling. 
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We have tested our disk-integration code by comparing 
its results with tho se of the code for horizontally homo- 
geneous planets by IStam et all ([2006D (see Appendix [A) ■ 
We have not compared it against other codes for mod- 
clling signals of hori z ontally inhomogeneou s planets (see 
e.g. [Ford et all 120011: lOaklev fc Cash! 120091; iTinetti et all 
2006), because these codes ignore polarisation, which is 
the most interesting feature of our code. A comparison be- 
tween calculated total fluxes would require us to disable 
the polarisation calculations in our adding-doubling code, 
since ignoring polarisation introduces er rors of up to sev- 
eral percent in total flux calculations (|Stam &; Hovenierl 
200$. From the comparison with the code for horizon- 
tally homogeneous planets applied to horizontally inhomo- 
geneous planets using the weighted sum approximation, 
in which weighted sums of flux vectors reflected by hor- 
izontally homogeneous planets are used to approximate 
the flux vectors of horizontally inhomogeneous planets (see 
Appendix|2l, we conclude that our code is accurate enough 
for application to horizontally inhomogeneous planets, pro- 
vided enough pixels are used across the disk, not only for 
resolving the spatial variations but also the variations in 
illumination and viewing angles across pixels. 

3. Sensitivity to horizontal inhomogeneities 

In this section, we present flux and polarisation (degree and 
angle) signals of planets with different types of horizontal 
inhomogeneities on their surfaces as calculated using our 
code to show the signatures of inhomogeneities and the dif- 
ferences with signatures of horizontally homogeneous plan- 
ets with similar surface albedos. Unless stated otherwise, 
each planet has a gaseous, Rayleigh scattering atmosphere 
with a total optical thickness of ~ 0.1 (no absorption), and 
a flat surface. For the surface albedo, we choose the ex- 
treme values of 0.0 (black) and 1.0 (white), since they give 
the largest contrast in flux and polarisation. 

Establishing which types of inhomogeneities could still 
be handled with e.g. a weighted sum approximation and 
which would need a full scale horizontally inhomogeneous 
approach with e.g. our code (from hereon: the Hi-code), is 
interesting for saving computing time when possible. The 
disadvantage of our Hi-code is namely the large amount of 
computing time it re quires as compared to the horizontally 
homogeneous code of IStam etTa l. (2006) (from hereon: the 
HH-code). For example, the Hi-code takes about 10 5 times 
more time than the HH-code for calculating the flux vectors 
of a horizontally homogeneous planet covered by a cloud 
layer with optical thickness 2, for phase angles a from 0° 
to 180° in steps of 2°, with the planet having been divided 
into pixels of 2° x 2°, and for a single wavelength. 

The difference in computation times would be irrele- 
vant if the computing time of the Hi-code were negligi- 
ble. Unfortunately this is not the case. This time depends 
strongly on the properties of the model atmosphere, espe- 
cially on the absorption and scattering optical thickness 
of the gases and particles in the atmosphere and on the 
angular variation of the single scattering properties of the 
scattering particles. In particular adding polarization to the 
flux calculations increases the time by at least an order of 
magnitude, and of course, the timing increases almost lin- 
early with the number of wavelengths at which calculations 
are required (since these calculations are independent of 
each other, they could easily be done in parallel). For the 



cloudy model planet calculations described in the previous 
paragraph, the HH-code takes about a minute on an aver- 
age workstation. 

The difference in computing time is not spent in the ra- 
diative transfer calculations themselves, since the two codes 
use the s ame adding-doubling radiative transfer algorithm 
(based on lde Haan et al1l!987f) and were run with the same 
numerical accuracy. Also, because both planets were hor- 
izontally homogeneous, the radiative transfer calculation 
had to be performed only once for each code. Indeed, in 
the Hl-code, the additional computing time is mostly spent 
in the integration of the flux vectors across the planetary 
disk. In particular, for a = 0°, the number of 2° x 2° pixels 
across the disk is more than 8000. For each pixel, the ap- 
propriate Fourier coefficients have to be determined from 
the list of calculated coefficients, they have to be summed 
to calculate the local reflection matrix, and with that the 
locally reflected flux vectors. Then, the locally reflected flux 
vectors have to be rotated to the planetary scattering plane 
in order to be summed up to calculate the reflected flux vec- 
tor of the planet. Even though the computing time per pixel 
can be relatively small, the mere number of pixels (which 
of course depends on the planetary phase angle) can result 
in long computing times. 

3.1. Planets with spots 

The first type of horizontally inhomogeneous planets have 
black surfaces with homogeneously distributed white spots. 
We start with spots consisting of single white pixels that 
cover 4% of the planet's surface, and let the area of each 
spot grow until the spots cover 68% of the surface. Figured] 
shows irF n and P s (see Eq. [3]) as calculated using the HI- 
code as functions of the planetary phase angle a for these 
planets. As can be seen, increasing the percentage of white 
pixels from 4% to 68%, increases irF n smoothly from 0.09 
to 0.58 at a = 0°, while P a decreases from 0.6 to 0.08 
at a = 90° (this phase angle does not coincide with the 
maximum of P s ). With increasing surface brightness, the 
peak value of P s shifts from a = 98° (4% white) to a = 116° 
(68% white). 

Figure [2] shows the difference between the results from 
the Hl-code and those from the HH-code (the latter com- 
bined with with the weighted sum approximation, as de- 
scribed in Appendix \K§. In Fig. [21 it is clear that when 
the planet is almost horizontally homogeneous (only 4% 
covered by white pixels), the maximum relative difference 
in flux, AirF n is very small (~ 0.64% for a = 0°), with 
the Hl-code giving the slightly higher values. With increas- 
ing coverage, AirF n increases, showing a maximum value of 
10% for 57% coverage in Fig. [21 to decrease again when the 
planet is almost homogeneously white. 

Because the degree of polarisation is itself a relative 
measure, we show the differences in P s between the Hl-code 
and the HH-code as absolute differences. As can be seen 
in Fig. [21 the absolute difference AP S is largest (~ —6% 
percent point around 90°) for the darkest planet and de- 
creases towards the whiter planets. The HH-code produces 
at (almost) all phase angles a larger P s than the Hl-code. 
Interestingly, the polarisation phase curve of the darkest 
planet is significantly less symmetrical when horizontally 
inhomogeneities are accurately taken into account (Hl- 
code) than with the HH-code and the weighted sum ap- 
proximation. 
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Fig. 1. irF n and P s as functions of a for black, cloud- 
free planets with white spots occupying 4% (red, dotted 
line), 21% (green, dashed line), 41% (blue, dashed-dotted 
line), 57% (grey, dashed-triple-dotted line) and 68% (pur- 
ple, long-dashed line) of the planet. 

3.2. Planets with center continents 

The second type of horizontally inhomogeneous planets 
have black surfaces and a circular "continent" of white pix- 
els at the center of the planetary disk facing the observer 
(rotation of the planet, hence the rotation of the continent 
in and out of the field-of-view is not taken into account 
here). Figure |3J shows irF n and P s as functions of a for 
planets with continents that cover from 5% to 98% of the 
disk. Both for the flux and the polarisation, the curves for 
80% and 98% coverage virtually overlap, because in these 
cases, the black pixels are all located along the limb of the 
planet and hardly contribute to the total signal. The flux 
phase functions in Fig. [3J have very similar shapes as those 
in Fig. [1] although the latter are darker for the same surface 
coverage of white pixels, which is not surprising since they 
have more black pixels on the front side of the disk. The 
polarisation phase curves in Fig. [3] clearly have different, 
more asymmetrical shapes than those in Fig. [TJ except for 
the largest coverages. The peak in the polarisation phase 
curves indicates the phase angle where the continent disap- 
pears into the planet's nightside. 

Comparing the flux and polarisation signals of the plan- 
ets with white continents as calculated with our HFcode to 



Fig. 2. Relative difference A7rF n and the absolute differ- 
ence AP S between the flux and polarisation phase func- 
tions of the white spotted planets of Fig. Q] as calculated 
using the HFcode and the HH-code. The coverage of white 
spots is 4% (red, dotted line), 21% (green, dashed line), 
41% (blue, dashed-dotted line), 57% (grey, dashed-tripple- 
dotted line) and 68% (purple, long-dashed line). 



the signals calculated using the HH-code (Fig. [4j), it is ob- 
vious that the differences A7rF n and AP S are smallest when 
the planetary disk is almost completely covered by the con- 
tinent. In Fig. |4l the largest difference in the flux is about 
33%, for a coverage of 32% and at a = 0°, with the HFcode 
giving the higher fluxes. The difference in polarisation, AP S , 
clearly shows the strong asymmetry around a — 90° of P s 
when calculated with the HFcode. 

Figure [5] shows A7rF n and AP S for the same planets 
except with a white surface and a black continent. Not sur- 
prisingly, for these planets, at most phase angles, TrF n is 
smaller when calculated with the HFcode than with the 
HH-code for each percentage of coverage at most phase an- 
gles, because of the concentration of black pixels in the 
centre of the disk. The difference AP S is less asymmet- 
ric than for the black planets with white continents (cf. 
Fig. [3]). The polarisation phase functions as calculated us- 
ing the two codes are thus similarly shaped. The maximum 
of the polarisation phase function, however, does depend 
strongly on the code, and is thus sensitive to the distri- 
bution of pixels across the disk, with the HFcode yielding 
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Fig. 3. irF n and P s as functions of a for black, cloud- 
free planets with a white continent on the center of the 
planetary disk facing the observer, for different coverages 
of the continent: 5% (black, solid line), 12% (red, dashed 
line), 21% (blue, dashed-dotted line), 32% (green, dashed- 
triple-dotted line), 80% (magenta, long-dashed line) and 
98% (gray, dotted line). 




planetary phase angle (deg) 

Fig. 4. AirF n and AP S between the phase functions of 
the white continent planets of Fig. [3] as calculated using 
the Hi-code and the HH-code. The coverage of the con- 
tinents is: 5% (black, solid line), 12% (red, dashed line), 
21% (blue, dashed-dotted line), 32% (green, dashed-triple- 
dotted line), 80% (magenta, long-dashed line) and 98% 
(gray, dotted line) of the planetary disk. 



much higher maximum values of P s than the HH-code, ex- 
cept for the smallest coverages. Indeed, for a white planet 
with a black continent covering 80% of the disk, P s is almost 
50% higher (around a — 80°) calculated with the Hi-code 
than with the HH-code. Note that for this coverage, the 
difference in flux is relatively small (~ —10%) and similar 
to that for a coverage of 5%. The difference in sensitivity 
to the spatial distribution of albedo across the planetary 
disk between flux and polarisation clearly illustrates the 
strengths of combined flux and polarisation measurements. 

3.3. Planets with hemispherical caps 

The third type of horizontally inhomogeneous planets have 
black surfaces and two, equally sized white caps on oppos- 
ing sides of the planet. We will present results for three 
locations of the caps: 0° (the caps cover the "north" and 
the "south" poles of the planet), 45°, and 90° (the caps 
are on the "eastern" and "western" sides of the planetary 
disk). Planets with caps at 45° are not mirror-symmetric 
with respect to the planetary scattering plane. Therefore, 



Stokes parameter U will usually not be zero. In the follow- 
ing, we will therefore use Eq.[2]to define the degree of linear 
polarisation instead of Eq. [3] 

First, we'll discuss the signals of planets with the caps 
covering the poles of the planets. Figure [6] shows TrF n and 
P for different coverages of the caps as functions of a as 
calculated using the Hi-code. The shape of the flux phase 
function is very smooth, and at a — 0° , irF n increases from 
0.07 for 23% coverage, to 0.56 for 77% coverage. In the 
latter case, the planet is basically white with a black, equa- 
torial belt (for a completely white planet, irF n would equal 
0.7, see Fig. [T5]). The polarisation phase function is quite 
symmetric for small caps, and P decreases from almost 0.55 
for 23% coverage to 0.08 for 77% coverage at a = 90° (note 
that the maxima of the polarisation phase functions occur 
at somewhat smaller and larger phase angles, respectively). 

In Fig. [7J we show A7rF n and AP for the planets in 
Fig.[S]when calculated using the Hi-code and the HH-code. 
The difference in flux is negative for all values of a: the 
Hi-code thus gives smaller fluxes than the HH-code for the 
same coverage. At a = 0°, AirF u is -0.14 for 23% cover- 
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Fig. 5. Similar to Fig. 01 except for white, cloud-free planets 
with black continents. The coverage of the continents is: 
5% (black, solid line), 21% (red, dashed line), 46% (blue, 
dashed-dotted line), and 80% (green, dashed-triple-dotted 
line) of the planetary disk. 



Fig. 6. irF n and P as functions of a for black, cloud- 
free planets with hemispherical caps occupying 23% (black, 
solid line), 43% (red, dashed line), 61% (blue, dashed- 
dotted line), and 77% (green, dashed-triple-dotted line) of 
the planet. The caps are located at the 'north' and 'south' 
poles of the planet. 



age, increases up to -0.21 for 43% coverage, and decreases 
again to -0.07 for 77% coverage. In differences in polarisa- 
tion show that for small coverages, the Hi-code gives sig- 
nificantly higher values of P than the HH-code. The reason 
is of course that for small coverages, the white pixels are 
mostly located near the limb of the planetary disk, where 
they contribute little to the total signal. 

Figure[S]shows 7rF n and P for planets with caps covering 
23% of the disk for three different locations of the caps. The 
three different flux phase functions show the relatively small 
effects of different fractions of the caps being visible at a 
given phase angle. The polarisation phase functions are all 
fairly symmetric around a = 90°, especially for the planet 
with its caps at the poles of the planet. 

Figure [9] shows A7rF n and AP for the three model plan- 
ets of Fig. [5] as calculated using the HH-code and the HI- 
code. As can be seen, the A7rF n are fairly independent of 
the caps' position angle. This implies that it would be dif- 
ficult to retrieve information about the position of the caps 
from the flux alone. Polarisation appears to be more sensi- 
tive to the location of the caps. In particular, the AP are 
several percentage points depending on the location and 
on a. Around a = 90°, where exoplanets have a relatively 



high chance of being observed with direct detection meth- 
ods, the differences in AP are about 0.12 for the planet 
with the eastern and western caps, implying that the ac- 
curacy of the polarimetry should be larger than 0.10 to be 
able to establish their existence. 



4. Application to Earth-like planets 

In this section, we present flux and polarisation signals of 
planets with horizontally inhomogeneous Earth-like surface 
coverages and patchy cloud layers. In reality, the Earth ex- 
hibits a large variation in atmospheric temperature and 
pressure profiles, cloud properties (both on macro- and 
micro-scales) and surface properties. To avoid introducing 
too many variables, we assume a single temperature and 
pressure profile across our mod el planet, the so-called m id- 
latitude summer profile from (McClatc hev et aT1ll972h . a 
single type of cloud particles and cloud properties (optical 
thickness and vertical distribution), and a surface that is 
covered by either ocean or sand (a continent). 
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planetary phase angle (deg) 

Fig. 7. A7rF n and AP S between the phase functions of the 
capped planets of Fig. \5\ as calculated using the Hi-code 
and the HH-code. The caps occupy 23% (black, solid line), 
43% (red, dashed line), 61% (blue, dashed-dotted line) and 
77% (green, dashed-triple-dotted line) of the planet. 

4.1. Surface in homogeneities 

First, we present the signals of a planet with a sandy conti- 
nent surrounded by ocean. The continent extends between 
longitudes of ±22° and latitudes of ±50° (measured with 
respect to the subobserver point). The sand surface reflects 
Lambertian (unpolarized and isotropically) with an albedo 
of 0.25. This albedo is taken from the ASTER spectral li- 
brary and should be representative for a sand surface on 
Earth at A = 0.55 /im. The ocean is black. The planetary 
atmosphere is cloud-free, with a total gas optical thickness 
of 0.1 (no absorption), which corresponds to a wavelength 
of 0.55 /xm. 

Figure [TU] shows irF n and P s as functions of a for the 
model planet as calculated with our Hi-code and the HH- 
code (combined with the weighted averages method), and, 
to compare, for model planets that are completely covered 
by sand or ocean. Both in ttF u and in P s , there are signif- 
icant differences between the signals from the Hi-code and 
the HH-code, depending on the phase angle. At a = 0°, the 
difference AttF s is ~ 30%. With increasing phase angle, the 
fraction of the continent on the planet's nightside increases, 
and irF n calculated with the Hi-code decreases faster than 
that calculated with the HH-code. The polarisation phase 
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Fig. 8. irF n and P as functions of a for black, cloud-free 
planets with white, hemispherical caps that occupy 23% of 
the planet. The locations of the caps are: 0° (black, solid 
line) , indicating the caps are on the 'north' and 'south' poles 
of the planet, 45° (red, dashed line), and 90° (blue, dashed- 
dotted line), with the caps on the 'eastern' and 'western' 
sides of the planetary disk. 



function calculated using the Hi-code is more asymmetric 
than that calculated using the HH-code. The largest abso- 
lute difference between the two curves is ~0.14 (14%) at 
a = 106°. 

Incidentally, for this planet, the flux and polarisation 
signals calculated using the Hi-code and the HH-code are 
virtually equal at a = 90°. At this phase angle, flux and 
polarisation observations would thus not help to establish 
the existence of the continent. However, interpreting flux 
observations at smaller (larger) phase angles using the HH- 
code would result in an overestimation (underestimation) of 
the coverage with sand, while interpreting only polarisation 
observations at smaller (larger) phase angles using the HH- 
code would result in an underestimation (overestimation) of 
the coverage with sand. The interpretation of the combined 
flux and polarisation observations using the HH-code would 
fail and thus reveal that assuming a homogeneous mixture 
of sand and ocean pixels is not realistic. 
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Fig. 9. A7rF n and AP S between the phase functions of the 
capped planets of Fig. [3] as calculated using the HFcode 
and the HH-code. The cap locations are: 0° (black, solid 
line), 45° (red, dashed line), and 90° (blue, dashed-dotted 
line) . 

4.2. Atmospheric inhomogeneities 

The next model planet has a black surface, a gaseous atmo- 
sphere with optical thickness of 0.1 (no absorption), and a 
patchy cloud layer with an optical thickness of 2.0. The 
clouds are composed of liquid water p articles described 
in siz e by the standard distribution of lHansen fc Travis! 
(|1974f ). with an effective radius of 2.0 fjja a nd an effective 
varian ce equal to 0.1 (model A particles of iKaralidi et al.l 
(2011)). The single scattering albedo of our particles at 0.55 
/j,m is 0.999534. For a detailed description of the single scat- 
tering properties of our cloud particles see IKaralidi et al.l 
(|201lh . The cloud layer is patchy (the clouds are composed 
of fully cloud covered pixels that cluster in random manner 
across the planetary surface) in the horizontal direction, 
but its vertical extension is the same all over the planet, 
namely from 3 km to 4 km in Earth's atmosphere (279 K 
to 273 K). 

Figure [IT] shows irF n and P as functions of phase an- 
gle a of the model planet for different cloud coverages as 
calculated using the HI- and the HH-code. The flux phase 
functions for the different cloud coverages have different 
strengths, but very similar shapes. The fase functions as 
calculated using the HH-code have slightly sharper fea- 



i.o : 

0.8 - 



ci- 




planetary phase angle (deg) 

Fig. 10. TrF n and P s as functions of a for a model planet 
covered by ocean and with a center continent of sand be- 
tween ±22° longitude and ±50° latitude, (green, dashed- 
triple-dotted line). Also plotted: the functions for a homo- 
geneous ocean planet (blue, dashed-dotted line), a homoge- 
neous sand planet (orange, solid line) and a weighted sum 
of these with ~16% sand and ~84% ocean (gray, dashed 
line). 



tures near a — 0°. Note that the bump in the curves near 
a = 30° is the signa ture of the primary rainbow (see e.g. 
IKaralidi et al.l 1201 ll) : light that is scattered once by the 
cloud particles. As expected, the difference between the 
fluxes calculated with the two codes are small for small 
and large percentages of cloud coverage. For an intermedi- 
ate, almost Earth-like, coverage of 42.3%, A7rF n is as large 
as 15% at a = 0°. In this case, using the HH-code to in- 
terpret the flux reflected by the patchy cloudy planet near 
a = 0° would yield a cloud coverage of ^53%. 

The polarisation phase functions clearly show the dif- 
ferent contributions of the light that is Rayleigh scattered 
by the gas molecules above and in particular between the 
patches of clouds (the strong values around a = 90°) and 
that of the light that is scattered by the cloud articles (the 
primary rainbow). Clearly, the larger the coverage of the 
clouds, the smaller the contribution of the purely Rayleigh 
scattered light to the total signal. Adding clouds to our 
model planet decreases the amount of purely Rayleigh scat- 
tered light, and thus decreases P around a = 90°. The 
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Fig. 11. irF n and P s as functions of a for a model ocean 
planet with patchy clouds that cover 16%, 42.3%, or 83.9% 
of the planet. The curves as calculated using the HH-code 
and the weighted sum approximation are also shown. 



strength of the primary rainbow in P is insensitive to the 
cloud coverage because adding clouds does not change the 
fraction of multiply scattered light within the clouds (the 
clouds have the same optical properties) , which mostly de- 
termines the strength of the rainbow on these planets. The 
Rayleigh scattering polarisation maximum around a = 90° 
does influence the contrast of the rainbow feature: for low 
cloud coverages it forms a shoulder on the Rayleigh scatter- 
ing maximum, while for high coverages, it is a local maxi- 
mum. 

For small and large percentages of cloud coverage, the 
differences in the polarisation phase functions due to using 
either the HI- or the HH-code are at most a few percent 
points around a = 90° (—4% for 16%, and -3% for 83.9% 
coverage). The differences are largest for the intermediate 
cloud coverage: about 10% for 42.3% coverage. In this case, 
using the HH-code to interpret the polarisation reflected by 
the patchy cloudy planet near a = 90°, would yield a cloud 
coverage of ^25%. The degree of polarisation calculated 
with the Hi-code, and thus AP too, depends not only on 
the cloud coverage, but also on the locations of the clouds 
across the planet. Not surprisingly, this sensitivity is highest 
for intermediate cloud coverages. 



Fig. 12. irF n and P s as functions of a for a model ocean 
planet with a central sandy continent and 42.3% cloud cov- 
erage as calculated using the Hi-code (red, dashed line) and 
the HH-code and the weighted sum approximation (black, 
solid line). For comparison, the curves for a model ocean 
planet without the continent and with 42.3% cloud cover- 
age are also shown (blue, dashed-dotted line). 



4.3. Atmospheric and surface in homogeneities 

Finally, we present the flux and polarisation signals of plan- 
ets with patchy clouds (see Sect. 14.21) . and a sandy conti- 
nent in the middle of a black ocean (see Sect. 14. ll) . Figure[T2l 
shows 7T-F n and P as functions of a for 42.3% cloud cover- 
age. The phase functions as calculated using the HH-code 
and the Hi-code results for a cloudy planet without a con- 
tinent are also shown. 

The continent strongly increases 7rF n . The increase will 
usually depend on the location of the clouds and the con- 
tinent. In this case, the continent is located in the middle 
of the planetary disk, and thus has a large influence. Since 
only a small fraction (~ 8%) of the pixels on the disk con- 
tain sand, the flux calculated with the HH-code is only 
marginally higher than that calculated with the HH-code 
and without a continent (see Fig. fTTj) . Inversely, if the HH- 
code would be used to interpret the flux signal of the cloudy 
planet with the continent, a 55% continental coverage and 
24% clouds would be found. 

The presence of the Lambertian reflecting continent be- 
low the clouds decreases the polarisation phase function and 
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shifts the maximum P towards larger a, when compared to 
the polarisation phase function of the cloudy ocean planet. 
The primary rainbow is still visible in the phase function, 
albeit more like a shoulder than a local maximum. 

Using polarimetry only, a straightforward fit to the po- 
larisation phase function with the HH-code would estimate 
the continental coverage at about 17% and the cloud cover- 
age at 24%. When both irF n and P are taken into account 
the best fit is acquired for the case of 22% continental cov- 
erage and 19% cloud coverage. 

5. Summary and conclusions 

We have presented a numerical code that can be used to cal- 
culate disk integrated total and polarised fluxes of starlight 
that is reflected by horizontally inhomogeneous exoplan- 
ets, e.g. planets that are covered by oceans and continents, 
and/or overlaid by a patchy cloud deck. 

For most types of model planets, our code for horizon- 
tally inhomogeneous planets (the Hi-code) is computation- 
ally much more expen sive than the code for horizontally ho- 
mogeneous planets bv lStam et al.l (|2006l) combined with the 
so-called weighted sum approximation to simulate fluxes 
of horizontally inhomogeneous planets (the HH-code). In 
the latter method, total and polarized flux signals of dif- 
ferent horizontally homogeneous planets are summed us- 
ing weighting factors depending on which fraction of the 
illuminated and visible part of the horizontally inhomoge- 
neous planet is represented by each type of planet. Only for 
planets with a large variation of horizontal inhomogeneities, 
such as a large number of different surface and cloud cover- 
ages, the computing time for the HH-code approaches that 
for the Hi-code. 

The main advantage of using the Hi-code instead of the 
HH-code is obviously the ability to simulate signals of hor- 
izontally inhomogeneous planets. Other advantages of the 
Hi-code are that since every pixel on the planet is treated 
separately, it allows taking into account effects of e.g. the 
non-sphericity of a planet, shadowing by planetary rings, 
and a spatial extension of the illuminating source. 

Since for most model planets, our Hi-cod e consumes 
much more computing time than the HH-code (|Stam et al.l 
2006) , it is interesting to investigate the influences of hori- 
zontal inhomogeneities on the total and polarized fluxes of 
starlight that is reflected by a planet. For three types of 
horizontally inhomogeneous planets covered by black and 
white surface pixels and overlaid by a gaseous atmosphere, 
we have calculated the fluxes and the degree of polarisa- 
tion as functions of the planetary phase angle a and com- 
pared them with results from the HH-code assuming the 
same percentage of black and white surface pixels and the 
same model atmosphere. Horizontal inhomogeneities can 
leave significant traces in both the reflected total flux irF n 
and the degree of linear polarisation. However, while hori- 
zontal inhomogeneities appear to mostly influence the total 
amount of reflected flux and not so much the shape of the 
planet's flux phase functions, they can strongly change the 
planet's polarisation phase functions in shape and strength. 
Indeed, fitting a planet's polarisation phase function using 
the HH-code would yield very different fractions of disk 
coverage (up to several tens of percent). 

We also used the Hi-code to calculate flux and polar- 
isation signals of Earth-like planets with surface and/or 
patchy clouds. These calculations confirmed that the shape 



of a planet's flux phase function is fairly independent of 
the horizontal inhomogeneities. Obviously, the absolute val- 
ues of the flux phase function does depend on the inhomo- 
geneities, but it will also depend on the radius of the planet 
(in this paper we assumed a planetary radius equal to one) , 
which thus would have to been known accurately to fit the 
flux phase function. The polarisation phase function ap- 
pears to be rather sensitive to the inhomogeneities, both 
its absolute values and the shape of the curve (e.g. loca- 
tion maximum value). Because the degree of polarisation 
is a ratio, it is independent of the radius of a planet. The 
flux and polarisa tion signals of planets are wavelength de- 
pendent (see e.g. lStamll2008t) . therefore the observability of 
horizontal inhomogeneities will also depend on the wave- 
length. In particular, at short wavelengths, the gas opti- 
cal thickness is larger than at longer wavelengths, and will 
thus hamper the observations of the surface. Future studies 
could focus on which spectral bands should be combined to 
optimize retrieval schemes. 

In the presence of liquid water clouds, the strength 
of the primary rainbow in the flux phase functions of 
Earth-like planets increases with increasing cloud coverage. 
Whether or not it would be detectable depends strongly on 
the sensitivity and stability of the observing instrument. In 
the polarisation phase function, the strength of the rainbow 
is almost independent of the cloud coverage as long as the 
planet's surface is very dark (i.e. does not add too much 
unpolarized flux to the total signal). A bright surface, such 
as a sandy continent, decreases the strength of the rainbow 
feature. 

Properly accounting for horizontal inhomogeneities ap- 
pears to significantly influence reflected fluxes and polarisa- 
tion signals, and should eventually be applied to interpret 
observations of horizontally inho mogeneous exoplan ets or 
e.g. observations of Earth-shine (|Sterzik et al.l 120121) . The 
HH-code (horizontally homogeneous planets combined with 
the weighted sum approximation), however, is still a strong 
tool for simulating signals to be used for the design and 
optimization of exoplanet observations, because its simula- 
tions cover the range of total and polarized fluxes that we 
can expect to observe. Because flux and polarisation phase 
functions have different sensitivities to the inhomogeneities, 
a combination of flux and polarisation observations would 
help to retrieve the actual planetary parameters. 

Our results, and especially the differences between the 
fluxes and degrees of polarization of the reflected starlight, 
indicate which accuracies should be reached with flux 
and/or polarization observations in order to detect hori- 
zontal inhomogeneities on a planet. As such they can drive 
the design of instruments for exoplanet characterization. 
Assuming a super-Earth exoplanet (with a radius of 1.5 r^) 
orbiting a Sun-like star at 1 pc from the observer at a 
40- meter telescope (such as the E-ELT), a back-of-the- 
envelope calculation shows that a few nights of integra- 
tion time (a total of 20 hours) could yield an accuracy of 
10~ 3 . The phase angle of a planet in an Earth-like orbit 
would change by ^3° during that time. For this calculation 
we ignored the influence of stellar background light on the 
observation, and we didn't include any actual instrument 
parameters, such as spectral bandwidths. Until results such 
as ours are combined with a realistic instrument and tele- 
scope simulator, the integration time estimate should thus 
be considered as a rough value. 
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Fig. 13. Normalised reflected fluxes irF n and irQ n as func- 
tions of the planetary phase angle a for a black planet 
with a gaseous atmosphere as calculated using the HH-code 
(solid lines) and the Hi-code (dashed lines). The solid and 
dashed lines coincide. 



A. Testing our numerical code 

Here, we present the results of testing our numerical code 
for horizontally inhomogeneous planets (from hereon the 
Hi-code). We have tested our code by comparing it to 
results of the cod e for horizontally homogeneous planets 
IStam et al.l (|2006l ) (from hereon the HH-code). 

A.l. Horizontally homogeneous planets 

For the first comparison between our Hi-code and the HH- 
code, we used both codes to calculate the light reflected 
by a model planet with a black surface and a gaseous, 
Rayleigh scattering atmosphere with a total optical thick- 
ness of about 0.1. For the Hi-code, the planet was divided 
into pixels of 2° x 2° (latitude x longitude). 

Figure [13] shows the excellent agreement between the 
reflected fluxes irF n and irQ n as functions of phase an- 
gle a calculated using the two codes (because the planet 
is mirror-symmetric with respect to the reference plane, 
7rC/ n = 0). The absolute difference between the fluxes cal- 
culated by the two codes is smaller than 10 -5 . Figure [TJ] 
shows the degree of polarisation P s (= —Q/F) (see Eq. [3]) 
calculated using both codes, and the absolute difference 
AP S between the curves. AP S is largest around small (15°) 
and large (160°) phase angles, which is due to the size of 
the pixels. Since in the Hi-code, we use one set of angles 9o, 
9, A</>, and /?2 for each pixel, the larger a pixel, the less this 
set of angles represents the range of angles across the pixel. 
The error is largest for pixels along the planetary limb and 
terminator, and thus for the large and small phase angles 
(at the small angles, the degree of polarisation in the centre 
of the disk is close to zero, making the contribution of the 
limb pixels thus significant). At a — 0° and 180°, P s of the 
planet is zero and thus AP S , too. 

For a model planet with a gaseous atmosphere, the dif- 
ferences between the results of the two codes appear to be 
fairly independent of the albedo of a (Lambertian reflect- 
ing) surface: with a white surface, the absolute differences in 
7rF n and 7rQ n , and AP S are similar to those with a black sur- 



Fig. 14. The degree of polarisation P s corresponding to the 
fluxes shown in Fig. [T3] (the two lines coincide) and the 
absolute difference between the two lines. 



face. For a model planet with a horizontally homogeneous 
cloud layer of optica l thickness 1, c o mpos ed of the model 
A cloud particles of iKaralidi et al.l (|2011l ). and a surface 
albedo of 0.1, the absolute difference in the fluxes increased 
slightly to a maximum value of 5xl0~ 5 and the maximum 
AP S increased to ~ 0.0062. These differences didn't change 
significantly when other (e.g. larger) cloud particles were 
used. 

The comparison with the HH-code shows that the inte- 
gration across the disk by our Hi-code is accurate enough 
for application to horizontally homogeneous planets, pro- 
vided small enough pixels are used across the disk. 

A, 2. Horizontally inhomogeneous planets 

For a second test of our Hi-code, we used model planets 
with surfaces covered by two types of 2° x 2° pixels that 
alternate in longitude and latitude (they thus look like 
chess boards). The equator of each planet coincides with 
the planetary scattering plane. The planetary atmospheres 
are gaseous, Rayleigh scattering and have an optical thick- 
ness of 0.1. While our Hi-code can handle t hese types o f 
planets, the HH-code cannot. Therefore, as in lStaml ([2008), 
weighted sums of flux vectors reflected by horizontally ho- 
mogeneous planets are used to approximate the flux vectors 
of horizontally inhomogeneous planets. The flux vector of 
a planet covered by J different types of pixels (with a dif- 
ferent atmosphere and/or surface) is thus calculated using 

J J 
irF(a) =^2w 3 nFj(a) with w 3 ■ = 1 (9) 

3=1 3=1 

with irFj the flux vector of starlight reflected by a planet 
that is completely covered by type j pixels, and with Wj the 
fraction of type j pixels on the horizontally inhomogeneous 
planet. 

Figure [15] shows 7rP n and P s as functions of phase angle 
a for a planet with its surface covered by alternating black 
and Lambertian reflecting white pixels as calculated with 
the Hi-code and the HH-code (the latter combined with the 
weighted sum approximation, see Eq. For comparison, 
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Fig. 15. irF n and P s as functions of a for a black planet 
(gray, dashcd-triple-dotted line), a white planet (blue, 
dashed-dotted line), and a planet with a gaseous atmo- 
sphere and its surface overed by alternating black and white 
pixels (black, solid line). The latter curves coincide with 
those calculated using horizontally homogeneous planets 
and a weighted sum approximation (red, dashed line). All 
planets are cloud-free. 



irF n and P s are also shown for horizontally homogeneous 
black and white planets. 

The absolute difference between the fluxes of starlight 
reflected by the black-and-white planet as calculated us- 
ing the Hi-code and the HH-code is smaller than 7xl0 -5 , 
and between the degrees of polarisation smaller than 0.005 
across the whole phase angle range. It appears that the 
differences between the results of the two codes depend 
slightly on the surface albedo: decreasing the albedo of pix- 
els from 1.0 to 0.24 (which is typical for a sand surface at 
0.55 Atm), the difference in the flux decreases to 5xl0 -5 
and in the polarisation to 0.004. 

The comparison with the HH-code applied to horizon- 
tally homogeneous planets and using the weighted sum ap- 
proximation shows that our Hi-code is accurate enough 
for application to horizontally inhomogeneous planets, pro- 
vided enough pixels are used across the disk. 
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